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Abstract 



Aims. This work aims to provide a new mass estimate for the Fornax cluster and the Fornax-Eridanus complex, avoiding 
methods like the virial or fits of X-ray emission profile, which assume that the system is in equilibrium, probably not 
the case of Fornax, still in process of formation. 

Methods. Our mass estimate is based on the determination of the zero-velocity surface which, in the context of the 
spherical infall model permits an evaluation of the total mass inside such a surface. The zero-velocity surface radius 
7?o was estimated either by a running median procedure or by fitting the data to the velocity field expected from the 
spherical model, including effects of the cosmological constant. The velocity field in a region within 20 Mpc of the 
Fornax center was mapped using a list of 109 galaxies whose distances have an average accuracy of 0.31 mag in their 
distance modulus. 

Results. Our analysis indicates that the mass of the Fornax cluster itself inside a radius of [2.62 — 5.18] Mpc is [0.40 — 
3.32] x 10 14 M while the mass inside [3.88 — 5.60] Mpc, corresponding to the Fornax-Eridanus complex is [1.30 — 
3.93] x 10 14 M Q . 

Key words, galaxies - clusters - individual: Fornax 
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In the hierarchical scenario of structure formation, clusters 
of galaxies are one of the largest structures observed in na- 
ture. Clusters have been assembled relatively late in the 
history of the universe, being located in the intersection of 
filaments constituting the cosmic web (Voit 2005). Clusters 
are bound essentially by the gravitational action of the so- 
called dark matter with the luminous (baryonic) component 
given only a minor contribution to the gravitational poten- 
tial of the system. Moreover, most of the baryonic matter 
is under the form of hot and warm gas filling the intra- 
cluster medium, which is detected by its X-ray emission. 
Clusters of galaxies have a particular interest in cosmology 
since the evolution of their number density above a given 
mass provide strong constraints on parameters characteriz- 
ing different world models (Rosati et al. 2002). 

Masses of clusters are generally estimated by using the 
virial relation, which presupposes that the system is dy- 
namically relaxed. This is not the case during phases in 
which the cluster undergoes a merger episode or accretes 
mass through filaments (see, for instance, Girardi & Biviano 
2002). Moreover, this method is affected by the eventual in- 
clusion of interlopers by projection effects in the considered 
sample and by dynamical friction, which may introduce an 
important bias between the velocity dispersion of galaxies 
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and dark matter particles. The reliability of the virial re- 
lation as a mass estimator has been checked by numerical 
simulations performed by different authors. Danese et al. 
(1981) and Fernley & Bhavsar (1984) concluded that pro- 
jection effects are important and may affect considerably 
the virial masses. Similar results were obtained by Perea 
et al. (1990), who reached the conclusion that the virial, 
in general, overestimates the masses unless interlopers are 
eliminated. A more recent study by Biviano et al. (2006) 
led to more quantitative results: the virial relation overes- 
timates true masses by about 10% if the simulated clusters 
have more than 60 members, with uncertainties increasing 
up to 50-60% for objects having 15-20 members. 

Besides the virial, other methods have been frequently 
employed as mass estimators of clusters. The X-ray emis- 
sion profile of the hot intracluster gas can be used to trace 
the gravitational potential of clusters, under the assump- 
tion of hydrostatic equilibrium (Sarazin 1986; Reiprich & 
Bohringer 2002). The accuracy of this approach was tested 
by cosmological simulations, which indicate that when 
masses are evaluated inside a radius at which the mean 
cluster density is about 500-2500 times the critical density, 
the uncertainties are of the order of 14-29% (Evrard et al. 
2006). The masses of galaxy clusters can be also estimated 
through the analysis of gravitational lensing, since the grav- 
itational field of clusters distorts the image of galaxies situ- 
ated behind them (Broadhurst et al. 1995). A comparison of 
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masses resulting from X-ray data and strong lensing shows 
that values drawn from the latter method are, on the aver- 
age, twice those obtained from the former procedure (Wu 
2000). These differences could be the consequence of an 
oversimplification of the lens model and/or the violation 
of the hypothesis of the gas isothermality. More recently, 
weak gravitational lensing, a technique permitting to track 
the gravitational potential of these objects by the distor- 
tion induced in the shape of background galaxies (see, for 
instance, Mellier 1999) has been used to estimate masses 
of different clusters. However, the presence of nearby fil- 
aments leads in general to an overestimate of masses by 
10-30% (Metzler et al. 1999). 

All the above mentioned methods correspond to scales 
less or of the order of 1-2 Mpc, which are typical dimensions 
of the central region of galaxy clusters. At larger scales, 
corresponding to the surroundings of clusters, galaxies are 
probably falling onto the cluster for the first time. These 
outskirt galaxies despite of being bound to the cluster are 
not in dynamical equilibrium. In this case, the knowledge 
of the velocity field of these objects may lead to an estimate 
of the central region mass. In fact, such an approach was 
proposed by Lynden-Bell (1981) and Sandage (1986) based 
on the spherical infall model. The motion of the outskirt 
galaxies is supposed to be essentially radial and from the 
knowledge of the distance Ro, at which the radial velocity 
with respect to the center of mass is zero, the mass inside 
such a surface can be estimated. The spherical model pre- 
dicts also the existence of caustics, surfaces at which (the- 
oretically) the galaxy number density is infinite (Rcgos & 
Geller 1989). The profile of the caustic amplitude, as seen 
in a phase-space diagram for the outskirt galaxies, can be 
used as a mass estimator with an accuracy of a factor of 
two (Diaferio 1999). 

In the present paper we intend to present a new esti- 
mate of the mass of the southern cluster located in Fornax 
(Abell S0373) at a distance of 20 Mpc. An early survey 
in the Fornax area performed by Ferguson (1989) indicates 
that probably 340 galaxies are cluster members and a fit of 
the projected density with a King profile suggests a core ra- 
dius of about 0.7°. Despite of being a cluster less rich than 
Virgo, the system presents different interesting features. 
Its main structure is centered in NGC 1399. According to 
Drinkwater et al. (2001), dwarf galaxies form a distinct 
population which is probably infalling onto the main sys- 
tem. Using the method by Diaferio (1999) mentioned above, 
which does't assume dynamical equilibrium, Drinkwater et 
al. (2001) estimated the projected mass inside a radius of 
1.4 Mpc as (7 ± 2) x 10 13 M Q . Since then a large amount 
of data on galaxies within the Local Universe, including 
radial velocities and precise distances have been accumu- 
lated. This justifies a novel study of the velocity field in the 
Fornax region as well as new estimates of its mass, based on 
the analysis of the velocity field of outskirt galaxies, mod- 
eled by the spherical infall model and avoiding problems 
present in other mass indicators as mentioned above. This 
paper is organized as follows: in Section 2 the available data 
is presented, in Section 3 mass estimates are discussed and 
finally in Section 4 the main results are summarized. 

2. The data 

In the past decade a significant number of galaxies present 
in the Local Volume had their distances measured with 



a quite good accuracy, in particular thanks to data ob- 
tained with the Hubble Space Telescope (Karachentsev et 
al. 2002a, 2002b, 2006, 2009). Besides galaxies present in 
the neighborhood of the Local Group a substantial effort 
was also made to increase the database on the Virgo cluster. 

Distances to galaxies in the Local Universe have been 
estimated from different methods: 

1) TRGB, based on the luminosity of the Tip of the 
Red Giant Branch, considered as one of the most effi- 
cient methods to determine distances of nearby galaxies, 
practically independent on their morphological type. The 
method requires images in two ore more photometric bands 
obtained with WFPC2 or ACS cameras on board of the 
HST, yielding an accuracy of about 7% on distances de- 
rived by such a procedure (Rizzi et al. 2007). A consol- 
idated list of distances for galaxies in the Local Volume 
is given in the Catalog of Neighboring Galaxies (hereafter 
CNG, Karachentsev et al. 2004). Galaxies from CNG with 
only TRGB or Cepheid distances were used, including some 
new determinations (Karachentsev et al. 2006, Tully et al. 
2006). 

2) The surface brightness fluctuation method (SBF), 
applied to early type galaxies, assumes that the old stellar 
population present in those objects gives the main contribu- 
tion to their luminosity. The method presupposes that the 
brightness distribution is not affected by irregularities as, 
for instance, that introduced by the presence of dust clouds. 
Using this approach, Tonry et al. (2001) determined SBF 
distances for 300 E and SO galaxies with typical errors of ~ 
12%. Galaxies in this sample are distributed over the whole 
sky, extending up to cz ~ 4000 km s^ 1 and having a median 
velocity of ^1800 km s _1 . 

3) Blakcslcc et al. (2009) undertook a two-color 
ACS/HST imaging survey including 43 early type galaxies 
situated in the Fornax core (the ACS Fornax Cluster Survey 
project, hereafter ACS-FCS), deriving SBF distances with 
errors of about 8%. To the ACS-FCS list were added 18 
dwarf ellipticals belonging to the cluster, having SBF dis- 
tances with an accuracy of about 9% estimated by Jerjen 
(2003) and Dunn & Jerjen (2006). These authors suggest 
from the S-shaped pattern distribution of these galaxies 
that Fornax is still in process of formation. 

4) Two galaxies within 15 Mpc of the Fornax cluster 
with distances measured with an accuracy of 5% by using 
SNIa light curves (Tonry et al. 2003) were also included in 
our database. 

5) Kashibadze (2008) determined distances for 402 edge- 
on spiral galaxies selected from the 2MASS Flat Galaxy 
Catalog (2MFGC, Mitronova et al. 2004), having radial ve- 
locities less than 3000 km s _1 . Using a multiparametric NIR 
Tully-Fisher relation, distances with an accuracy of about 
20% were obtained. The zero point of the luminosity-line 
width relation was established by using 15 galaxies with 
distances derived from cepheids and TRGB data. 

6) We supplemented the aforementioned samples with 
a compilation of distances by Tully et al. (2008, 2009) de- 
rived from the Tully-Fisher relation calibrated for optical 
(B, V, R, I) magnitudes. This compilation relays on numer- 
ous HI line and photometric observations carried out by 
Methewson & Ford (1996), Haynes et al. (1999), Tully 
& Pierce (2000), Koribalski et al. (2004), Springob et al. 
(2005), Theureau et al. (2006) and other authors. Again, 
the zero point of these relations were set by using a sam- 
ple of 40 galaxies with distances determined by cepheids 
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Figure 1. Sky distribution in equatorial coordinates of 562 galaxies within 20 Mpc from the Fornax center. Markers of 
different diameters indicate three distinct distance ranges from the Local Group. The circle in the lower panel corresponds 
to the zero velocity sphere projection (Rq = 4.60 Mpc) for the Fornax-Eridanus complex. 



and TRGB data. Finally, the list of galaxies by Springob et 
al. (2007) (SFI++ sample), not considered by Tully et al. 
(2009) was also used in our compilation. 

The first four samples are referred as precise data be- 
cause their typical errors don't exceed 10-12%, while the 
last two datasets are mentioned as Tully-Fisher data. 

In this paper we will generally follow Karachentsev & 
Nasonova (2010) and examine the 3D sample considering 
galaxies with limited spatial distances from the Fornax clus- 
ter center. As it was discussed by the authors, this approach 
isn't free of systematical effects because galaxy distances 
are measured with errors and their significance is different 



at the proximate and the distant boundary of the spherical 
volume (the so-called Malmquist bias; see Figure 2). 

Our initial list includes 1140 galaxies within 30 Mpc 
of the Fornax cluster (available in the electronic version of 
the paper). However, in the present study, we will be fo- 
cused in a region of 20 Mpc around the center of the clus- 
ter, representing a sample of 562 objects which, in prin- 
ciple, would permit an estimate not only of zero velocity 
surface but also of the transition region between bound 
and unbound objects, the latter tracking essentially the 
Hubble flow. The characteristics of these galaxies are given 
in Table 1. The first column indicates the distance esti- 
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Table 1. Database for galaxies within 20 Mpc from the 
Fornax center. 



Distance method 


Cm 


N 


G 


TRGB + Ceph 


0.15 


113 


7.1 


SBF (Tonry) 


0.25 


69 


3.3 


ACSFCS+ 


0.19 


55 


3.9 


SNIa (Tonry) 


0.10 


2 


1.4 


TF (IR) 


0.40 


69 


2.1 


TF (opt) 


0.40 


254 


4.0 


All 




562 


27.8 



mate method; the second column gives the mean error in 
the distance modulus; the third column gives the number of 
galaxies and the last column gives the sample goodness de- 
fined as G = (AVIOO) 1 / 2 -^ 1 (Kudryaet al. 2003). Figure 1 
shows the projected distribution of these galaxies in the sky. 
Galaxies are marked as circles and their diameters indicate 
different distance ranges. 

2.1. The velocity field 

In order to map the velocity field around the Fornax center 
and to have an estimate of the velocity-distance relation, ac- 
curate radial velocities and distances are required. Then, in 
the next step, these data must be converted into distances 
and velocities relative to the cluster center. 

Data on radial velocities, mostly from HI observations, 
were obtained in most cases from the same sources of dis- 
tances. When not available, NED data on heliocentric ve- 
locities were used. Observational errors in radial velocities 
are quite small (1-2 km s -1 ) in the case of HI observations 
(Tifft & Huchtmeier 1990) and can be neglected in compar- 
ison with distance errors (AV <C AR ■ Ho) in the scales of 
the nearest clusters (R > 15 Mpc). 

The transformation of heliocentric velocities into the 
Local Group reference frame was performed with the stan- 
dard apex parameters (Karachentsev & Makarov 1996) 
adopted in NED. If ip is the angular separation between 
the apex and a galaxy, then the converted velocity is 
V C on — V+V apex cos ip and the error of this conversion is not 

more than [(AV) 2 + (AV apex ) 2 + {Aip ■ V apex ) 2 ] 1/2 , where 
AV apex — 5 km s _1 and ip w 1%, so Aip- V apex w 3 km s _1 . 
Thus, the errors introduced by this transformation are 
about 6 km s _1 , being negligible in comparison with dis- 
tance errors. 

The gravitational effect of the Fornax cluster can be 
seen directly from the radial velocity vs. distance relation 
as an S-shaped wave. Radial velocities and distances rel- 
ative to the Local Group centroid for 98 galaxies in the 
cluster core (6 < 5°) are represented in the top panel of 
Figure 2. Here, precise distances for most galaxies were 
obtained within the special survey ACS-FCS with HST 
(Blakeslee et al. 2009). The centroid of galaxies forming 
the "virial column" at [19.58 ± 1.25] Mpc is marked by 
gray (see Section 3.1 for discussion of the Fornax clus- 
ter barycenter position). The plotted value of virial ra- 
dius, ±1.25 Mpc, is an approximate estimate based on 
the i?o value for the Fornax-Eridanus complex obtained 
in this paper (4.60 Mpc). The S-shaped curves correspond 
to a Hubble flow perturbed by a point-like masses of 
1.30 • 10 14 M Q and 3.93 • 10 14 M Q as the limiting cases within 



the confidence range (see Section 3.2 for details) in the case 
of line-of-sight passing exactly through the cluster center. 
The typical distance error bars for datasets (2), (3) and (4) 
are shown. 

The distribution of radial velocities and distances for 
remaining galaxies of the sample in periphery of the Fornax 
cluster (5° < 9 < 30°) is shown in the bottom panel of 
Figure 2. Here, the S-shaped lines having lower amplitudes 
describe the behavior of perturbed Hubble flow at angular 
distance 9 = 5°. The typical distance error bars for datasets 
(1), (6) and (7) are presented. 

However, a serious source of uncertainty is caused by 
the absence of data on tangential velocities. In order to 
estimate galaxy velocities with respect to the center of mass 
of Fornax, some model should be used and thus, the results 
turn out to be model-dependent. Taking into account the 
lack of data on the true velocity vector of galaxies, there are 
at least two approaches to obtain such a transformation. 

The first one assumes that the prevailing motion, which 
involves most of the galaxies under study is the asymptotic 
Hubble relationship (the model of the minor attractor). The 
second approach considers that galaxies are within the in- 
fall zone (the model of the major extended attractor), i.e., 
they do not follow the Hubble flow but instead are falling 
towards the cluster center. Both cases were discussed in 
details by Karachentsev & Nasonova (2010); see Figure 3 
sketching the relative positions of the considered galaxy 
(G), the observer (LG) and the cluster center). 

When a galaxy is located strictly in front or behind the 
cluster center (i.e. the angles A and 9 are small), both ap- 
proaches yield the same infall velocity toward the cluster 
center. When A is close to 90°, in the second case Vi n — > oo 
leading to a significant discrepancy between the two ap- 
proaches. However, there are no dramatic differences be- 
tween both methods in the Hubble diagram. Some galaxies 
move along the vertical axis appreciably, but the behavior 
of running medians (see next section) traces the infall of 
galaxies towards the cluster in a similar way. Nevertheless, 
as we shall see later, the second method yields systemati- 
cally slightly larger values of Rq. The scatter of galaxies in 
the Hubble diagram also increases in the second case. These 
considerations suggest the first approach to be preferred. 

Finally, in order to reduce the role of the unknown tan- 
gential component of the velocity and to avoid further un- 
certainties, we decided to select for our analysis only galax- 
ies situated approximately in front and behind the cluster, 
i.e. in a cone satisfying the conditions A < 45° or A > 135°. 
The number of galaxies satisfying this additional constraint 
within 20 Mpc of the Fornax center is 164 and their pro- 
jected distribution in the sky is shown in Figure 4. 

It should be mentioned that the role of possible chaotic 
tangential velocities of the galaxies had been studied by 
Karachentsev & Kashibadze (2006). They have performed 
numerical simulations, adding some random tangential 
component to the observed radial velocity. Their model- 
ing of the Hubble flow in the vicinity of the Local Group 
showed that typical tangential velocities with amplitudes 
of 35 and 70 km s _1 produce a statistical uncertainty in 
the evaluation of the zero-velocity surface radius as small 
as ±2% and ±4% respectively. 
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Figure 2. Radial velocity vs. distance relation for galaxies in the Fornax cluster region with respect to the Local Group 
centroid. Galaxy samples with distances derived by different methods are marked by different symbols. The inclined 
straight line traces Hubble relation with the global Hubble parameter Hq — 73 km s _1 Mpc -1 . The virial zone is 
filled with gray. Two S-shaped lines correspond to a Hubble flow perturbed by masses of 1.30 • 10 14 M Q (dotted) and 
3.93 • 10 14 M Q (dashed) as the limiting cases within the confidence range. The typical distance error bars for each dataset 
are shown. Top: the cluster core within angular distance 8 < 5°, the S-shaped lines indicate the expected infall at 9 = 0°. 
Bottom: peripheric regions with 5° < 6 < 30°, the S-shaped lines indicats the infall at 8 — 5°. 




Figure 3. Relative positions between the 
observer LG, the considered galaxy G and 
the cluster center C. Left panel: Case of a 
pure Hubble flow. Right panel : Case of a 
pure infall towards the Fornax center. 
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Figure 4. Spatial distribution of galax- 
ies around the Fornax cluster. The zones 
where A < 45° or A > 135° are blue shaded 
whereas the exclusion zones are white. 
The same color code indicates galaxies in- 
cluded or not in the analysis. For a ques- 
tion of clearness, lines of equal A for values 
in the range 15°- 90° are also shown. 



3. Mass estimates 

The difficulties with the different mass estimators were al- 
ready mentioned. In particular, the presence of interlopers 
lead to an overestimate of the mass ranging from 10% up 
to 60% when the virial relation is used, whereas estimates 
based on X-ray data have uncertainties ranging from 14- 
30%. 

The case of the Fornax cluster is quite particular. Unlike 
clusters as Coma or Virgo, which have a well defined center 
and a more or less smooth mass distribution, despite the 
presence of some substructures, the Fornax cluster has a 
complicated mass distribution, including many substruc- 
tures. Such a complexity is sketched in Figure 5 where 
it is possible to identify: the Fornax cluster centered on 
NGC 1399, the Eridanus subcluster, including the groups 
NGC 1407, NGC 1332, NGC 1395 and NGC 1398 as well as 
other small neighboring groups. It should be emphasized 
that the virial radii of these groups do not overlap, suggest- 
ing that these objects are separated and constitute gravita- 
tionally bound structures. However these groups overlap in 
the scale of the zero- velocity radius, indicating that despite 
the absence of dynamical equilibrium, they are probably 
bound to the main structure. 



3.1. The spatial position of the Fornax-Eridanus bar/center 

In the present paper we regard the dynamical center of the 
Fornax cluster (NGC 1399 group) situated near the core 
of the hot X-ray emitting gas (Jones et al. 1997, Scharf 
et al. 2005, Machacek et al. 2005) as the gravity center 
of the whole Fornax-Eridanus complex. Its spatial posi- 
tion was calculated as the mean position of all NGC 1399 
group members, yielding D = 19.6 Mpc, a = 3 h 32 m 30 s and 
S = — 35°51'15". For the moment this choice (as the first 
approximation) seems to be reasonable since the NGC 1399 
region has been investigated most detaily. The next step is 
to determine the spatial position of the Fornax-Eridanus 
complex barycenter as the weighted mean for position vec- 
tors of all the virialized substructures forming the complex. 
This will be possible after obtaining more observational 
data in a wider area of the complex. Still, the same calcula- 
tion techniques are applied by us for both center positions 
(with values R = 21.1 Mpc, a = 3 h 33 m 58 s , S = -28°44'45" 



adopted for the second case), and the resulting R values 
didn't differ significantly. The barycenter position associ- 
ated with NGC 1399 group yields 0.23 Mpc lower i?o value 
from precise data and 0.19 Mpc higher value from Tully- 
Fisher data, leading to the upper bound of discrepancy 
~ 0.2 Mpc. Generally speaking, the zero-velocity surface 
method is rather stable in the sense of a barycenter posi- 
tion (Karachentsev & Kashibadze 2006). 

Table 2 represents groups, triplets and pairs form- 
ing the Fornax-Eridanus complex with individual distance 
estimates (Karachentsev & Makarov 2008; Makarov & 
Karachentsev 2009, 2011). Columns of the table contain: 
(1) name of the brightest galaxy in a group/triplet/pair; (2) 
equatorial coordinates of the brightest galaxy for triplets 
and pairs or mean equatorial coordinates for groups (at the 
J2000.0 equinox); (3) number of galaxies in a system; (4) 
number of galaxies with measured distances; (5, 6) mean 
velocity with respect to the Local Group centroid and its 
error; (7) radial velocity dispersion; (8) harmonical radius; 

(9) integrated luminosity of the group in the K s band; 

(10) virial mass of the group; (11, 12) mean distance of 
the group and its error; (13) Hubble distance calculated as 
D h = H V hG where H Q = 73 km s" 1 Mpc^ 1 ; (14) and (16) 
angular distance to the Fornax-Eridanus complex barycen- 
ter for the cases of (14) dynamical center of the NGC 1399 
group situated near the center of the hot X-ray emitting 
gas and (16) dynamical center of all substructures forming 
the complex; (15) and (17) angle between the directions 
to the Local Group and to the Fornax-Eridanus complex 
barycenter for both cases mentioned above, respectively. 

The possible role of barycenter position is illustrated by 
Figure 6 representing velocity vs. distance diagrams where 
only the centers of the groups forming the Fornax-Eridanus 
complex are shown. The top panel represents radial ve- 
locities and distances of the groups relative to the Local 
Group centroid. The distance to the Fornax-Eridanus com- 
plex, 21.1 Mpc, corresponds to the spatial position of the 
complex barycenter as the weighted mean for position vec- 
tors of all the virialized substructures forming the complex. 
The sets of dashed and dotted lines indicate the Hubble 
flow perturbed by a point-like mass of 2.16 • 10 14 M© and 
0.87- 10 14 M Q respectively which corresponds to the mass of 
the whole Fornax-Eridanus complex estimated in this pa- 
per and the virial mass estimation of the Fornax cluster 
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Table 2. Groups, triplets and pairs around the Fornax-Eridanus complex with individual distances. 



Group 


RA, Dec 


N 


N D 


v hG 


± 


a v 


Rh 




lgM 


D 


± 


D h 


e 


A 


& 


A 




hhmmss ddmmss 






km/ 


s 


km/s 


kpc 


L 


Mq 


Mpc 


Mpc 






O 


o 


(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


NGC 1097 


024621.0- 


314430 


4 


1 


1189 


81 


140 


43 


11.26 


12.53 


15.9 


3.2 


16.3 


11 


131 


10 


142 


NGC1187 


030238.6- 


230819 


4 


1 


1315 


28 


48 


106 


10.60 


11.90 


18.4 


3.7 


18.0 


15 


96 


9 


125 


NGC 1201 


030408.0- 


275549 


2 


2 


1633 


48 


48 


312 


10.95 


11.94 


21.0 


0.8 


22.4 


12 


65 


7 


88 


NGC 1232 


030945.5- 


212514 


3 


2 


1582 


49 


70 


475 


11.07 


12.93 


19.2 


0.5 


21.7 


16 


86 


10 


113 


NGC 1291 


031718.6- 


425331 


3 


3 


703 


70 


99 


226 


11.03 


13.09 


8.2 


0.5 


9.6 


6 


170 


13 


159 


NGC 1302 


031723.8- 
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Figure 5. Fornax cluster (NGC 1399 group), Eridanus 
subcluster (NGC 1395, NGC 1407, NGC 1332, NGC 1398 
groups) and other neighboring groups of galaxies forming 
the Fornax complex (Makarov & Karachentsev, 2011). In 
the scales of virial radii (dark-blue circles) these overden- 
sities constitute separate bound groups, but they overlap 
substantially in scales of Ro (light-blue regions). The ve- 
locities of groups Vlg are indicated near the names of their 
central galaxies. NGC 1399 as the most massive galaxy in 
the complex situated near the center of the hot X-ray emit- 
ting gas is indicated as a purple circle marked with a cross. 
Black dots represent galaxies with measured distances. The 
large bold circle of radius 13.2° outlines the zero-velocity 
sphere (Ro = 4.60 Mpc) around the Fornax-Eridanus clus- 
ter bary center. 



itself (Makarov & Karachentsev 2011). Both sets include 
S-shaped lines signing different angular separations from 
the complex center (0°, 5°, 15° and 30°). The error bars 
indicate rms uncertainties in distances and velocities. 

The velocity vs. distance diagram relative to the cen- 
ter of the Fornax-Eridanus complex is shown in the bottom 
panel of the Figure 6. Two types of markers, circles and 
squares, correspond to the center of the hot X-ray emitting 
gas and the dynamical center of all substructures form- 
ing the complex, respectively. The solid line indicates the 
Hubble flow perturbed by a point-like mass of 0.87T0 14 M Q 
(the virial mass estimation of the Fornax cluster itself) 
while the dashed and dotted lines are constructed as least 
square estimates for both adopted center positions men- 
tioned above (0.53-10 14 M Q and 1.66-10 14 M© respectively). 

The structural complexity of the galaxy distribution 
suggests that probably the best mass indicator for this clus- 
ter should be based on the velocity field of the outskirt 
galaxies. Here two approaches will be adopted. In the first, 
the radius of the zero velocity surface Ro will be estimated 
from a running median procedure. From the knowledge of 
Ro, the mass inside such a surface follows immediately. In 
the second, the expected radial velocity profile derived from 
the spherical model is fitted to the data and again, the mass 
inside the zero-velocity surface results from the derived fit 
parameters. 

3.2. The running median 

The zero-velocity radius Ro can be estimated from a run- 
ning median procedure, using directly observational data 
(Figure 7). 

In order to account errors in observed distances and 
velocities and to give some estimates of the uncertain- 
ties, a Monte- Carlo simulations technique was used. For 
any galaxy i in a given dataset j, a corresponding dis- 
tance from the center is generated according to the relation 
Ri = Ri i0 bs + ARiftij, where Ri t0 bs is the observed distance 
from the Fornax center, ARi is an observational error as- 
sociated to the distance and /Jy is a normally distributed 
random number (with a = 1). Generated velocities were de- 
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Figure 6. Top: Radial velocities and distances of the groups relative to the Local Group centroid. The distance to the 
Fornax-Eridanus complex, 21.1 Mpc, corresponds to the spatial position of the dynamical center of all substructures 
forming the complex. The sets of dashed and dotted lines indicate the Hubble flow perturbed by a point-like mass of 
2.16 • 10 14 M Q and 0.87 • 10 14 M Q respectively, signing different angular separations from the complex center (0°, 5°, 15° 
and 30°). The error bars indicate rms uncertainties in distances and velocities. Bottom: The velocity vs. distance diagram 
relative to the center of the Fornax-Eridanus complex. Round and square markers correspond to the center of the hot 
X-ray emitting gas and the dynamical center of all substructures forming the complex, respectively. The dotted, solid 
and dashed lines correspond to the Hubble flow perturbed by a point-like mass of 0.53 • 1O 14 M , 0.87 • 1O 14 M , and 
1.66 • 1O 14 M respectively. 



rived by a similar procedure. About 10,000 datasets were 
generated and for each of them the running median method 
was applied, yielding different values of i?o, which were 
then averaged. The distribution of 10,000 i?o are shown 
in Figure 8 as well as the mean and median values and the 
90 %error band (for the window w = 1 Mpc) . The result- 



ing median values and corresponding errors for (A) minor 
attractor and (B) major extended attractor cases are given 
in Table 3. The hrst two columns correspond to the cases 
when only precise or Tully-Fisher data were used for Monte- 
Carlo simulations while the third column corresponds to 
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Figure 7. Hubble flow in the Fornax cluster reference frame for 109 galaxies with i?pc < 20 Mpc. The black solid 
polygon curve traces the running median on observational data with a window of 1 Mpc. The dotted and dashed curves 
correspond to a Hubble flow perturbed by masses of 1.30 • 10 14 M Q and 3.93 • 10 14 M Q respectively as the limiting cases 
within the confidence range. Only the galaxies with A < 45° or A > 135° are presented. Top: the case of almost pure 
Hubble flow ( "minor attractor" approach) . Blue curves trace the running medians on simulated data with a window of 1 
Mpc and the corresponding 90% error band. Bottom: the case of almost pure infall towards the Fornax cluster ("major 
attractor" approach). Coloured curves trace the running medians on simulated data with different window values. 



the whole dataset. The lines of the table correspond to the 
different median windows. 

With the Fornax cluster distance of ~ 20 Mpc and the 
typical uncertainty of 10 % for the ACSFCS+ galaxies pop- 
ulating the central core of the cluster, the observational dis- 
tance errors in the virialized zone will be of about 2 Mpc. 
Assuming the virial radius of the Fornax cluster to be ~ 1.5 



Mpc we should conclude that the cluster core galaxies can 
possibly distort the pattern in the Rq region for the Fornax 
cluster itself, but not for the whole Fornax-Eridanus com- 
plex. Anyway, the galaxies in the Fornax cluster core, i.e. 
those with Rpc < 3.5 — 3.7 Mpc, were not regarded in the 
Monte-Carlo data simulations, resulting to 109 galaxies (see 
Figure 7). 
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Figure 8. Distribution of 10,000 7? re- 
alizations simulated using Monte-Carlo 
technique, their mean and median values 
and the 90 % error band. The window used 
for the running median procedure is w = 
1 Mpc. 



Table 3. Radius Rq and its 90% error band (in Mpc) obtained from Monte-Carlo simulations. 



Window 
(Mpc) 


Precise distances 
(samples 1-4) 


TF distances 
(samples 5-6) 


All 

samples 


(A) (B) 


(A) (B) 


(A) (B) 


0.8 


4.48 4.56 
3.88-5.46 3.94-5.63 


4.87 5.15 
3.95-6.18 4.08-6.68 


4.41 4.67 
3.83-5.39 3.93-5.76 


1.0 


4.58 4.72 
3.87-5.50 3.95-5.75 


5.04 5.39 
3.95-6.32 4.09-6.81 


4.42 4.78 
3.82-5.42 3.93-5.79 


1.2 


4.51 4.66 
3.88-5.56 3.92-5.85 


5.15 5.51 
3.95-6.33 4.16-6.86 


4.39 4.83 
3.82-5.48 3.93-5.88 



The substantial dip in the running median of simula- 
tions based on precise data can be explained partially by the 
enormous velocity of NGC 1400 (see, for example, Perrett 
et al. 1997). Belonging to the NGC 1407 group, NGC 1400 
has an extremely low velocity relative to the Local Group 
(558 km s- 1 ) while NGC 1407 has V LG = 1771 km s" 1 . 
Both galaxies have roughly the same distance (26.4 Mpc for 
NGC 1400 and 28.8 Mpc for NGC 1407 according to Tonry 
et al. 2001). In the velocity-distance diagrams NGC 1400 
has a significant negative velocity relatively to the Fornax- 
Eridanus complex center, —380 km s _1 and —988 km s _1 
for the cases of minor and major attractors respectively, 
and therefore it appears in the legend region at the Fornax- 
centric distance of 9.6 Mpc shifting the running median 
downwards. 

According to the Table 3, the median estimate of i?o 
is 4.60 Mpc with confidence interval corresponding to 90 % 
error band of [3.88 - 5.60] Mpc. 

Once Rq is known, the mass inside the zero- velocity sur- 
face can be computed. Using the spherical model, including 
the effects of the cosmological constant, the mass inside i? 
is (E. Shaya, private communication; Karachentsev et al. 
2007) 

where 

'^>-i^-5(i^-"^is:- 1) - ,2) 

with n m being the mass density parameter. 

Figure 9 shows the ratio of masses of a galaxy system 
computed in the ACDM model and in the empty Universe 
model respectively as a function of f2 m . As one can see, 
the adopted uncertainty in A value affects as about 3% 
in mass estimation, which is negligible as compared with 
uncertainties caused by observational errors. 

ACDM parameters can be determined from 
WMAP data with a sufficient accuracy, i.e., 
H = 73.2±f% km s" 1 Mpc" 1 and O m = 0.238 ±0.015 




Figure 9. Ratio of masses of a galaxy system computed in 
the ACDM model and in the empty Universe model respec- 
tively as a function of Q m . The value of Q m = 0.238 ±0.015 
adopted in the ACDM model is marked with a vertical gray 
strip. 

(Spergel et al. 2007). Substituting these values into the 
Equation (1), we get 

(M T /M )o. 24 - 2.23 • 10 12 (i? /Mpc) 3 . (3) 

Notice that from the numerical solution of the equations 
describing the Lemaitre-Tolman model, modified in order 
to include effects of the cosmological constant, Peirani & 
de Freitas Pacheco (2008) obtained for the mass inside the 
zero- velocity surface 

M = 4.24 x 10 12 h 2 Rl M Q (4) 

The numerical coefficient was obtained for Ct m = 0.3 and 
for h — 0.73 both relations derived either analytically 
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Table 4. Masses resulting from the fit of eq.[5]to the present 
data. 



Parameter 


"Minor attractor" 


"Major attractor' 


Mass (in 10 i4 M ) 


0.93 


1.55 




0.38-2.10 


0.43-4.54 


Ro (in Mpc) 


3.44 


4.08 




2.57-4.52 


2.67-5.85 


a (in km s _1 ) 


300 


345 


a c (in km s _1 ) 


84 


190 



or numerically agree quite well. For Rq = 4.60 Mpc the 
corresponding total mass of the Fornax complex is M = 
2.16 x 10 14 M Q and for the limiting values of 90 % error band 
[3.88 - 5.60] Mpc we obtain M = [1.30 - 3.93] x 10 14 M Q 
as a confidence interval. 

3.3. Mass estimate from the radial velocity profile 

In this section, the mass inside the zero- velocity surface is 
estimated by fitting a theoretical profile directly to data. 
We follow the procedure developed by Peirani & de Freitas 
Pacheco (2006, 2008), who have numerically computed the 
velocity field of outskirt galaxies, based on the spherical col- 
lapse model including effects of the cosmological constant. 
This approach assumes that: a) most of the cluster mass 
is concentrated in the core, in which the shell crossing has 
already occurred and b) that orbits of galaxies outside the 
core are mainly radial. From the numerical calculations by 
Peirani & de Freitas Pacheco (2008), the velocity-distance 
relation is 

nfl)= , 3 ™ oR _^(™y" +,V3 . (5) 

Here M is the core cluster mass, R is the distance of 
the member galaxy to the cluster center, V(R) is the radial 
velocity of the galaxy with respect to the mass center, Hq is 
the present value of the Hubble parameter and n — 0.627. 
The relation above results from a fit of numerical data and 
is valid for the present time [z — 0), since it varies as a 
function of the redshift. Notice that eq. H]can be deduced 
from this equation when V(Rq) = 0. 

In principle, if accurate radial velocities and distances 
are available, from the fitting of eq. [3] to data it is possi- 
ble to derive both the core cluster mass and the Hubble 
parameter (Peirani & de Freitas Pacheco 2008). Here face 
to the uncertainties in the determination of Fornax-centric 
velocities, as explained in the previous section, we adopted 
the procedure of fixing Hq as 73 km s _1 Mpc -1 and com- 
puting the mass either for the "minor attractor" or for the 
"major extended attractor" models. 

Masses resulting from our fitting procedure are given in 
Table 4 for the two models adopted to estimate Fornax- 
centric velocities. The zero velocity surface radius Rq and 
the velocity dispersion a with respect to the mean infall flow 
are also given. The second column corresponds to velocities 
estimated from the "minor attractor" approach whereas the 
third column gives values derived from the "major extended 
attractor" model. The corresponding confidence intervals 
for mass and Rq values are estimated as 90 % error band. 
The third and the fourth lines give respectively the uncor- 
rected and the corrected velocity dispersions with respect 



to the infall flow towards Fornax. Notice that this method 
leads to values of masses that are, on the average, 50% 
smaller than the "running median" approach, being rather 
a mass estimate for the Fornax cluster core than for the to- 
tal Fornax-Eridanus complex. It is worth mentioning that 
the derived 1-D velocity dispersion with respect the bulk 
motion, (300-345 km s~ 4 ), is one order of magnitude higher 
than that observed in the vicinities of the Local Group but 
it drops to (84-190 km s _1 ) when the distance measure- 
ment errors are taken into account. Figure 10 shows the 
best fit between the expected velocity profile (eq. [5]) derived 
from the spherical model and data, whose velocities were 
computed according to the "minor attractor" model (upper 
panel) or to the "major attractor" model (lower panel). The 
number of galaxies used in the fitting procedure is limited 
by i?max satisfying the relation R msK H = V(R msK ) and 
actually accounts for 22 and 23 galaxies respectively. 

Considering the mass values estimated from these pro- 
cedures we adopt for the Fornax cluster itself a mass of 
1.24 x 10 14 M Q with the confidence interval of [0.40 - 
3.32] x 10 14 M©, and for the total Fornax-Eridanus com- 
plex a mass of 2.16 x 10 14 M & with the confidence interval 
of [1.30 — 3.93] x 1O 14 M0, which corresponds to a value 
half order of magnitude the Virgo cluster mass. In Table 5 
we compare previous mass estimates for the Fornax cluster 
(normalized to Dp or = 20 Mpc and Dsn = 25 Mpc) with 
the present study, based on the velocity field of outskirt 
galaxies. 

4. Concluding remarks 

The distribution of galaxies in the vicinity of the Fornax 
cluster indicates a substantial degree of subclustering that 
forming the Fornax complex. As a consequence, mass esti- 
mates based on the virial relation are affected by the fact 
that the system has not yet reached an equilibrium state 
and is still probably in formation (Dunn & Jerjen 2006). 
In fact dwarf galaxies have distinct dynamic properties in 
comparison with bright galaxies and are probably infalling 
into the system (Drinkwater et al. 2001). These difficul- 
ties may be circumvent by studying the velocity field of 
outskirt galaxies, which permits an estimate of the zero- 
velocity surface and, consequently of the mass inside such 
a surface as proposed originally by Lynden-Bell (1981) and 
Sandage (1986). 

In this investigation, a culled sample of 109 galaxies 
with measured distances and within 20 Mpc from the clus- 
ter center were used to study the velocity field in the neigh- 
borhood of the Fornax cluster. Since tangential velocities 
are unknown, in order to estimated the galaxy velocities 
with respect to the cluster center two main assumptions 
were made: in the first, it was supposed that the velocity 
vector is essentially dominated by the Hubble flow ( "minor 
attractor" model) while in the second, it was assumed that 
galaxies are infalling radially ("major attractor" model). 
The zero- velocity radius was derived by two different meth- 
ods: the "running median" and by fitting directly the data 
to the expected velocity profile derived from the spherical 
model, including the effects of a cosmological constant. 

The best fit radius of the zero-velocity sphere for the 
Fornax-Eridanus complex is estimated by us to be i?o 
4.00 Mpc with the confidence interval of [3.38 — 5.60] Mpc 
while the mass inside such a surface is M to t = [1.30 — 3.93] x 
10 14 Mq. At the distance to Fornax cluster 20 Mpc, the ra- 
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Figure 10. Velocity field in the Fornax region, for 164 galaxies with distances less than 20 Mpc from the cluster centre. 
The 90 % error band resulting from Monte-Carlo simulations is shown. Upper panel: velocities derived according to the 
"minor attractor" approach. Lower panel: velocities derived according to the "major attractor" model. 



dius Rq = 4.60 Mpc corresponds to 13.2° shown in Figure 5. 
Notice that within this circle there are almost all systems 
identified by Makarov & Karachentsev (2011) as virialized 
groups bound to Fornax cluster. Their total virial mass, 
M = 1.92 x 10 14 Mq, agrees with M tot , meaning that prob- 
ably only a small part of the Fornax- Eridanus complex mass 
is spreaded between the groups. 
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